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Abstract 

Conformal embedding of closed-universe models in a de Sitter back- 
ground suggests a quantisation condition on the available conformal time. 
This condition implies that the universe is closed at no greater than the 
10% level. When a massive scalar field is introduced to drive an infla- 
tionary phase this figure is reduced to closure at nearer the 1% level. In 
order to enforce the constraint on the available conformal time we need to 
consider conditions in the universe before the onset of inflation. A formal 
f | ' series around the initial singularity is constructed, which rests on a pair 

of dimensionless, scale-invariant parameters. For physically-acceptable 
models we find that both parameters are of order unity, so no fine tuning 
is required, except in the mass of the scalar field. For typical values of 
the input parameters we predict the observed values of the cosmological 
^ ' parameters, including the magnitude of the cosmological constant. The 

OS , model produces a very good fit to the most recent CMBR data, predicting 

C — . a low-^ fall-off in the CMB power spectrum consistent with that observed 

^ ' by WMAP. 

t-H ■ 

^ ■ 1 Introduction 

' Experimental evidence is now strongly in favour of the idea that a non-zero 

cosmological constant, or some form of 'dark energy', is currently responsible for 
around 70% of the total energy density of the universe |J . There have been many 
, theoretical attempts to justify the reintroduction of a A-term, something that 

was once viewed as a rather ugly and unnecessary extension to classical general 
relativity. Many of these are based on ideas from particle physics, including 
concepts such as spontaneous symmetry breaking and false vacua. These ideas 
Jv>( , can successfully explain the size of the cosmological 'constant' required during 

■ an inflationary phase, but all such models have great difficulty in explaining 

the current scale of the cosmological constant without implausible levels of fine 
tuning. 

It could be argued that the current scale of the cosmological constant sug- 
gests that its role may be more geometric, rather than field-theoretic. This idea 
fits in well with the gauge-theoretic viewpoint of gravity that we have developed 
elsewhere IB]- In the absence of any matter, a non-zero cosmological constant 
implies that the universe should be described by a de Sitter space. This space 
should then form the background for the gauge theory treatment of gravity, 
though the details of such a theory remain to be fully worked out. 

Here we explore some simpler consequences of the view that A should be 
viewed as a genuine constant with a geometric role to play. The ideas outlined 
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here are discussed in more detail in [3]. We first introduce a new picture of 
de Sitter geometry, which generalises the Poincare disk picture of hyperbolic 
geometry E| • Since all cosmological models are conformally equivalent (they 
all have vanishing Wcyl curvature), all cosmologies with a de Sitter end state 
will have a standard embedding in the conformal picture. These embcdding- 
s are revealing, and suggest a 'preferred' cosmology — one with closed spatial 
sections, in which the total available conformal time is ir/2. Such a model places 
the initial singularity symmetrically in the centre of the de Sitter picture. The 
requirement on the total conformal time suggests the operation of a quantisation 
condition, and we speculate as to how such a condition may arise. Of course, 
there are many theoretical reasons for preferring spatial flatness which is, after 
all, supposed to be one of the main predictions of inflation. For example, there 
are difficulties in constructing a homogeneous stress-energy tensor for Dirac 
fermions in anything other than a spatially-flat universe Qj- But closed 
universe models have their own attractive features, particularly in their inherent 
finiteness. Furthermore, in closed (or open) cosmologies the scale factor can be 
determined absolutely from the equation 



the scale factor through reheating in order to compare physical scales today 
with those during inflation. 

Several other authors have of course been interested recently in closed uni- 
verse models, motivated partly by the apparent \ow-£ turn down in the WMAP 
results 1 j. For example, Efstathiou 7 proposed a phcnomological model in 
which the low-fc primordial spectrum in a closed model has an exponential cut- 
off on a scale comparable to the curvature scale. On the other hand, Starobin- 
sky [S] , predicted that the effect of closure for adiabatic scalar perturbations was 
in general to increase low-^ values. This is clearly an area where more detailed 
calculations are necessary. 

On the observational side, Efstathiou (e.g. 0) and others have queried 
whether foreground effects on power spectrum estimation could be affecting the 
low-£ modes. However, it is likely that the low value seen for the quadrupole 
mode is still statistically significant, and thus represents something which needs 
explaining in the standard model. Other explanations which have been given 
include ^Oj which has a period of fast-roll in a flat model, and also topological 
attempts at explanation, via effects in a compact closed universe 

What we will look at here is simple closed universe, with no topological 
effects. But, by pursuing our boundary condition, we will find strong constraints 
among cosmological parameters. In particular, for a cosmological model with a 
given equation of state, the condition that the total elapsed conformal time is 
equal to ir/2 singles out a unique trajectory in the (CIm> ^a) plane. For example, 
a simple dust model predicts a universe that is closed at around the 10% level. 
This represents an upper limit, and both inflationary and radiation-dominated 
epochs drive this figure down. A prediction of the cosmological model can 
equally be interpreted as a prediction of A and, for typical input parameters, 
we find that 
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This result simplifies calculations for perturbations, as we do not have to track 



A « exp(-6N)l. 
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Here N is the number of e-foldings in the inflationary region, and turns out to 
be roughly 46, putting A in precisely the observed range. 

The model of inflation we consider is the simplest available — that of a 
massive scalar field. This produces an extremely tight model that agrees well 
with all current experimental data. In order to apply our boundary condition we 
need to study the evolution of the scalar field from the initial singularity, through 
the inflationary region, before matching onto a standard cosmological model. It 
has become quite common in the inflationary community to ignore the initial 
singularity and concentrate instead on suitable initial conditions for inflation. It 
might have been hoped that the peculiar nature of the inflationary stress-energy 
tensor, which violates the weak energy condition, might circumvent some of the 
standard singularity theorems. But this is not the case. Inflationary models are 
geodesically incomplete in the past, and the cosmological equations can be easily 
run backwards in time to reveal the initial singularity |12l I13| . In this respect, 
inflation has little to say about the 'specialness' of the initial singularity (zero 
Weyl curvature) which can be regarded as one of the outstanding problems in 
physics |14| . 

Expanding the field equations around the initial singularity can be performed 
using an iterative scheme in powers of t and Int. The resulting series are gov- 
erned by two parameters, which effectively control the degree of inflation and 
the curvature. Fixing both of these parameters to be of order unity produces 
inflationary models in a closed universe which are consistent with observation. 
This appears to provide a counterexample to statements that it is difficult to 
obtain closed universe inflation without excessive fine tuning |151 1161 1171 IT5| . 

As the universe exits the inflationary region it evolves as if it had started 
from an effective big-bang, with a displaced time coordinate. Photons have 
travelled an appreciable distance by the end of inflation, which alters how we 
apply the w/2 boundary condition. The result of these effects is the imposition 
of a see-saw mechanism linking the current state of the universe and the initial 
conditions. The more we increase the number of e-folds during inflation, the 
smaller the value of the cosmological constant, and vice- versa. With initial 
conditions chosen to give the required number of e-foldings to generate the 
observed perturbation spectrum, we find that the predicted universe is closed 
at the level of a few percent, in excellent agreement with observation. More 
detailed calculation also reveals a dip in the low-£ part of the power spectrum. 
However, these calculations are difficult and we comment on the problems that 
must be overcome in order to make this prediction robust. 

Unless stated otherwise we work in units where G = c = h = 1. Where it 
adds clarity, factors of G are included, so that G has dimensions of (distance) 2 . 

2 Conformal pictures of de Sitter space 

De Sitter space is a space of constant negative curvature, forming the Lorentzian 
analogue of the non-Euclidean geometry discovered by Lobachevskii and Bolyai |SJ 
[3]. Two-dimensional non-Euclidean geometry has an elegant construction in 
terms of the Poincare disc (which was, in fact, first given by Beltrami in 1868, 
and later rediscovered by Poincare In the Poincare disk model, geodesies 

are represented as 'd-lines' — circles that intersect the disc boundary at right- 
angles (see figure ^1 . Here we develop a similar picture for 2-dimensional de 
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Figure 1: The Poincare disc. Points inside the disc represent points in a 2- 
dimensional non-Euclidean (hyperbolic) space. A set of d-lines are also shown. 
These are (Euclidean) circles that intersect the unit circle at right angles. Given 
a (Mine and point A not on the line, one can find an infinite number of lines 
through A that do not intersect the line. 



Sitter space, which is sufficient to capture the key features of the geometry. We 
start with an embedding picture, representing de Sitter space as the 2-surface 
defined by 

T 2 - X 2 - Y 2 = -a 2 , (3) 

where (T, X, Y) denote coordinates in a space of signature (1,2) and a is a 
constant. The resulting surface is illustrated in figure [21 The figure illustrates 
a key property of the entire de Sitter geometry, which is that spatial sections 
are closed, whereas the timelike direction is open. So de Sitter space describes 
a closed universe that lasts for infinite time. One can set up local coordinate 
patches for which spatial sections are flat or open, but these coordinates are not 
global. These are discussed in the following section. Null geodesies are straight 
lines formed by the intersection of the surface and a vertical plane a distance a 
from the timelike axis. Despite the fact that the space is spatially closed, the 
furthest a photon can travel is half of the way round the universe. 

The preceding picture is the 'standard' one of de Sitter space, though it is not 
the one originally put forward by de Sitter |19l I2(J| . Further spaces of constant 
curvature can be created by topological identifications, and among these is the 
antipodal map on spatial sections. In terms of embedding coordinates (T, X 1 ), 
this is the map (T, X 1 ) i— > (T, —X 1 ). The result has spatial sections of the form 
IP™ as opposed to S n , where MP" denotes the real n-dimensional projective 
plane. The projective plane can be quite hard to visualise. For example, RP 2 
involves taking a Mobius strip, and applying a further twisted identification on 
the 'long' side (figure 01. The resulting manifold in non-orientable, and cannot 
be embedded in three dimensions without self-intersection (it is similar to a 
Klein bottle). Of more physical relevance is KP 3 , which is orientable, though 
the space contains non-contractable loops and is equally hard to visualise. MP 3 
is also a group manifold, so as a background space it does have some attractive 
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Figure 2: Two-dimensional de Sitter Space. The timelike direction is vertical, 
and spatial sections are closed. The right-hand diagram shows a null geodesic, 
which is a straight line in the embedding space. 



(a) (b) 

Figure 3: (a) A representation of the projective plane. Common arrows denote 
surfaces to be identified, together with the orientation, (b) The Mobius strip, 
shown for comparison. 
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future infinity 




past infinity 



Figure 4: Conformal diagram of de Sitter space. The vertical axis represents 
conformal time, running from to n. This version of the diagram is appropriate 
for a 4-dimcnsional spacetime, with 8 and <fi coordinates suppressed. 



properties. The 3-sphere S 3 is the group manifold of the three dimensional 
spin group (and also SU(2)), and M.P 3 is correspondingly the group manifold 
of SO (3). From our perspective, the standard framework of S n spatial sections 
is the one of greatest interest, though this is largely an aesthetic judgement. 
There does not appear to be any obvious physical criteria for favouring one over 
the other. 

There are a range of conformal representations one can adopt for de Sitter 
space |21|. Probably the simplest is the Carter-Penrose diagram, shown in 
figure In this picture the temporal coordinate represents conformal time, 
and the spatial section represents the effective radial coordinate. A 2-sphere is 
suppressed at each point. If one wants to represent the simpler 2-dimensional 
de Sitter space, the diagram needs to be wrapped onto a cylinder. While this 
picture is instructive, we seek a picture that embodies some of the features 
of the Poincare disk model. So, for example, geodesies should have a simple 
representation in terms of curves of constant (Minkowski) distance from some 
point. To achieve such a picture, we start by considering the spatial section at 
T = 0. This section is a ring of radius a, which is mapped onto a straight line in 
a (1, 1) Lorentzian space via a stereographic projection. Null geodesies from this 
section are represented as 45° straight lines in Lorentzian space. Since geodesies 
from opposite points on the ring meet at infinity, we arrive at a boundary in the 
timelike direction defined by a hyperbola. This construction provides us with 
a Lorentzian view of de Sitter geometry. Timelike geodesies in de Sitter space 
are represented by hyperbolae that intersect the boundary at a (Lorentzian) 
right-angle (see figure [5J. 

There are many fascinating geometric structures associated with de Sitter 
geometry, many of which mirror those of non-Euclidean geometry For example, 
one can always find a reflection that takes any point to the origin [221 • One can 
then prove a number of results at the origin, where the geodesies are all straight 
lines, and the results are guaranteed to hold at all points. A similar approach 
can be applied to anti-de Sitter space, except now the diagram is rotated through 
90°, as it is the timelike direction that is formed from a stereographic unwrap- 
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Figure 5: The Lorentzian view of two-dimensional de Sitter space. The bound- 
ary is defined by two hyperbolae (shown with thick lines). All geodesies through 
the origin are straight lines, and null geodesies are always straight lines at 45°. 
Two further geodesies, one spacelike and one timelike, are also shown. These 
are hyperbolae which do not pass through the origin. The timelike geodesic in- 
tersects the boundary in such a way that the two tangent vectors have vanishing 
Lorentzian inner product. 



ping of a circle. A valuable mathematical tool for studying these geometries is 
provided by conformal geometric algebra pfl I2H| . which is essentially a Clifford 
algebra for the conformal space of two dimensions higher. So a 4-dimensional dc 
Sitter spacetime is studied in the Clifford algebra of a 6-dimensional space with 
signature (2,4). Points, lines (including geodesies), planes and spheres are all 
represented by graded multivectors in the algebra, and union and intersection 
operations can be defined analogously to projective geometry [221 - 

The metric associated with the Lorentzian view of 2-dimensional de Sitter 
geometry has the conformally-flat structure 

ds 2 = — 7 % — (dt 2 -dx 2 ). (4) 

(a 2 +x 2 -t 2 ) 2V ' v ; 

This makes it clear that null geodesies must remain as straight lines in the x-t 
plane. A similar line element is appropriate for four dimensional spacetimes, but 
these are not commonplace in the literature (except for the simple case of flat 
cosmologies) . We finish this section by outlining the steps required to transform 
a standard FRW line element in the form 

ds 2 = dt 2 - + ^ 2 2 /4)2 (dr 2 + r 2 (d6 2 + sin 2 (0) dcf> 2 )) (5) 

into the spacetime-conformal line element 

ds 2 = — (dt 2 - dr 2 - x 2 {de 2 + sin 2 (9) d(j> 2 )), (6) 

where / is a (dimensionless) function of t and r. 
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The first problem to address is that the r coordinate in equation (JSJ is 
assumed to be dimensionless. To rectify this we introduce a constant distance 
A and replace the line clement of JSJ) with 

ds " = df2 - (tf+kr'y ( dr * + r 2 (d6 2 + sin 2 (0)dtf)). (7) 

Here t, r, A and R all have units of distance (assuming c — 1). As usual, the 
constant k is either ±1 or zero. Comparing the angular terms in equations (JJJ 
and (jnjl we see immediately that 



(8) 



r _ 2XrR(t) 
/ ~~ A 2 + kr 2 ' 

The coordinate transformation from (t, r) to (t, r) must satisfy 

r 

dt = f cosh(w) dt-\ — sinh(u) dr 
r 
r 

dr = f sinh(u) dt H — cosh(it) dr. (9) 
r 

For flat cosmologies (k — 0) we simply set t/A equal to the conformal time rj, 
where 

dt' 



*=lwr (10) 

and we then have r = r and the hyperbolic angle u is set to zero. For non-flat 
cosmologies it is perhaps surprising to find that u is non-zero. That is, there is a 
mismatch between the conformal coordinate frame and the cosmological frame. 
It follows that the conformal time rj is not the same as the time-like conformal 
coordinate t. 

Using the integrability conditions for the coordinate transformations, and 
setting the initial singularity to t = 0, the following solutions are found in the 
three cosmological scenarios 0]: 

f = —7 — r = q \ —t, — — rr ;] ~ closed 

J Rsin{rj) J \X 2 + r 2 -t 2 J A 

/ = ^ flat (11) 

t 1 -( 2At \ t 

f= RsJnW)= g U 2 +t 2 -r 2 ) A ^ 

Here g and g are functions (in general elliptic) satisfying differential equations 
depending on the matter content. 



3 A boundary condition for conformal time 

All cosmological models are conformally flat, and can all be interchanged via 
conformal transformations. Furthermore, all models containing a cosmological 
constant, and which do not recollapse, will tend towards a de Sitter end state. 
Such models should fit neatly within the conformal diagram of figure [S] with 
future infinity represented by the upper hyperbolic boundary. A flat section 
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Figure 6: Flat and open cosmological models, de Sitter space contains sections 
representing both flat and open spacetimes. A flat spacetime (left) consists 
of the space contained inside a light-cone located on past infinity. An open 
spacetime fits inside the lightcone from the origin. Both pictures can be used 
to illustrate cosmological models, with the initial singularity represented by a 
hyperbolic spacelike surface. An example of these surfaces is shown on each 
diagram. 



of de Sitter space corresponds to a region contained within a light cone from a 
point located on the boundary at past infinity (see figure^ . Surfaces of constant 
cosmic time are represented as hyperbolae, any one of which can then be chosen 
to represent the initial singularity. Similarly, an open section of de Sitter space 
is represented by the area inside a light-cone from a point in the middle of the 
de Sitter picture. An open A-cosmology will have the initial singularity located 
on a spacelike hyperbola. 

The diagrams for flat and open cosmologies make it clear that one is not 
employing the full de Sitter geometry in a symmetric manner. Closed models, 
on the other hand, can be given a more natural embedding, since their initial 
singularity can be represented by a spacelike surface placed anywhere in the 
conformal diagram of figure El For all but one choice of position of the initial 
singularity, the future asymptotic de Sitter state of the model will not match 
onto the future infinity boundary of the conformal diagram. However, suppose 
that we insist that t = represents the initial singularity. In this case the entire 
history of the universe can be conformally mapped into the top half of the de 
Sitter diagram. We suggest that this gives a natural boundary condition, as the 
centre of the diagram is the most natural, symmetric place to locate the big- 
bang singularity. (This choice also hints at a possible pre-big bang phase, but 
this is not explored here.) It might be thought that our proposal just amounts 
to a choice of coordinates. However, we now demonstrate that it has physically 
measurable consequences, and also that it can be recovered via an eigenvalue- 
type condition on the underlying equations, hinting at its origin as a possible 
quantisation condition. 

One immediate physical consequence of our proposal is that any photon 
emitted from the initial singularity may travel a maximum of precisely one- 
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quarter of the way round the universe in the entire future evolution of the 
universe. An alternative way of saying this is that the past horizon projected 
back to t = should cover half of the de Sitter geometry. If we let </> denote 
an equatorial angle on a 3-sphere, a photon travelling round the equator will 
satisfy 

1-5. < 12 > 

where t is cosmic time. Traversing the entire universe corresponds to running 
from < 4> < 2n. If a photon is only to travel one-quarter of the way round we 
therefore require that 

oa dt_r ca <m 

~R~ J + AiP/3 + SirpIP/3) 1 /* ~ ' ( ' 

That is, the total available conformal time is ir/2. This constraint can also be 
arrived at through an alternative route that works entirely within the conformal 
representation of cosmological models [32]. A de Sitter space centred on t = 
= t has the conformal line element 

12A 2 

ds 2 = — — — (dt 2 - dx 2 - dy 2 - dz 2 ) . (14) 

A(A 2 + r 2 — t 2 ) 2 v v ' 

Clearly the only forms for / in equation Ijlll) that have any chance of matching 
onto this final state are those for a closed universe. Furthermore, the function 
g must satisfy 

A 2 A\ 1/2 1 



where 



But since 



lim g(x) = HH - 

X^oo V 3 / X 

2 At . . 

X = A 2 + r 2 _ t 2 = tan W' ( 16 ) 

.9 = (17) 
Hsm{rj) 

we must then have Rcos(r]) tending to a constant at large times. This is only 
possible if r\ tends to tt/2, recovering our earlier boundary condition. This 
derivation is instructive in that it reveals how the constraint can be imposed as 
a straightforward boundary condition on a differential equation. In this case, 
the equation for g(x) is 

2/1 , 2n f d 9\ 2 . d 2 8irG\ 2 p A 2 A 

X(1 + X) UJ + dx- {gx) = ^^ + — ■ (18) 

The task then is to solve this equation subject to the boundary condition that 
g falls off as 1/x for large x- Viewed this way the constraint can be thought of 
as a 'quantisation condition' applied as the universe is formed, which one might 
hope would emerge from a quantum theory of gravity. 

In order to understand the implications of our boundary condition, we turn 
to considering flow lines in the (Um,^a) plane (see figure 0). The plots show 
a shows a series of flow lines starting from Qm = 1, = 0, refocusing around 
the spatially flat case, VLm + = !• One can show that, for a large range 
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of initial conditions, by the time we reach the current value of Am ~ 0.3 most 
models are not far off spatial flatness. This prediction contrasts with models 
without a cosmological constant, where any slight deviation from the critical 
density in the early universe is scaled enormously by the time we reach the 
present epoch, implying that the parameters in the early universe are highly 
fine-tuned. The presence of a cosmological constant therefore goes some way to 
solving the flatness problem on its own, without the need to invoke inflation. 

Each of the flow lines in figured has a different value for the total evolved 
conformal time. Imposing a constraint on this picks out a single preferred 
trajectory, resulting in the two curves shown in figure |H1 (one for matter-filled 
and one for radiation- filled). As the universe is expected to be matter dominated 
for most of its history, the solid line in figure El is the more physically relevant 
one. Taking the present-day energy density to be around Am = 0.3, we see 
that fl\ = is predicted to be around Q\ « 0.83. That is, a universe that is 
closed at around the 10% ratio. Such a prediction is reasonably close to the 
observed value, though it is ruled out by the most recent experiments [Q. In 
order to improve the prediction, we need to use up a greater fraction of the 
conformal time before we enter the matter-dominated phase. Such a process is 
also necessary to solve the horizon problem, and the simplest means of achieving 
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u 0.2 0.4 0.6 0.8 1 

Q M 

Figure 8: Critical paths as predicted by the de Sitter embedding. The solid line 
represents a matter-dominated universe, and the broken line shows radiation, 
for comparison. The straight line is the critical case of spatial flatness. By the 
time we reach Qm ~ 0.3 the universe has been driven close to critical density. 



this is via an inflationary phase. 

4 Scalar field inflation 

One can argue about the extent to which inflation really does solve the original 
problems it was intended to. For example, we have seen that the presence of 
a cosmological constant alone goes a long way to solving the flatness problem, 
without requiring inflation to put the universe in an effective de Sitter phase. 
But when it comes to generating structure on scales greater than the horizon, 
inflation is currently the only option. The simplest model of inflation models 
the matter content as a real, time-dependent, homogeneous massive scalar field. 
For a range of initial conditions this system shows the expected inflationary 
behaviour. The evolution equations for this model are 

A + H2 ~l + ^(20 2 --V) = O (19) 

and 

4> + 3#0 + m 2 (j> = 0. (20) 

Given a solution to the pair of equations (|19Jl and l|20|l . a new solution set is 
generated by scaling with a constant <r and defining 

H'(t) = aH(at), <j>'(t) = cj>{ot), m! = am, A' = a 2 A. (21) 

This scaling property is valuable for numerical work, as a range of situations 
can be covered by a single numerical integration. A range of physical quantities 
are invariant under changes in scale, including the conformal time r\ as can be 
seen from equation l|l(J|) . The scaling property does not survive quantisation, 
however, so has to be employed carefully when considering vacuum fluctuations. 
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The initial conditions for equations H19|l and (|2U|I are usually set at the start 
of the inflationary period, where they are viewed as arising from some form of 
quantum gravity interaction. But in order to apply our boundary condition we 
need to track the equations right back to the initial singularity, as this is the 
only way that we can keep track of the total conformal time that has elapsed. 
Evolving the inflationary equation backwards in time is entirely justified, as 
we do not expect quantum gravity to play a role until much earlier in the 
history of the universe. Inflation on its own does not eliminate singularities from 
cosmology |12l 113 ) . Furthermore, by specifying initial conditions around the 
singularity, the states of the fields at the onset of inflation are fully determined 
by a pair of parameters, making the model highly predictive. 

As the universe emerges from the big bang the dominant behaviour of H is 
to go as 1/ (3t) . Equation 119fl then implies that </> must contain a term going 
as ln(i). But this in turn implies that H must also contain a term in tln(t), 
in order to satisfy equation H20J) . Working in this manner we conclude that a 
series expansion in powers of t and ln(t) is required to describe the behaviour 
around the singularity. At this point it is convenient to define the dimensionless 
variables 

t rn 
u=-, M=— (22) 
ip Trip 

with tp and m p the Planck time and mass respectively. The series expansion 
about the singularity at t = can now be written 

1 oo 1 CO 

l P i=0 l p i=0 

which ensures that the expansion coefficients are all dimensionless. Substituting 
these into the two evolution equations, and setting each coefficient of ln(u) to 
zero, we establish that 

dHo 2 uA 8iru ( dcf>o \ Wircf) 

Hi = —u— uH Q ' 



du 3 3 V du ) 3 du 

3u 3 



(24) 



with further algebraic equations holding for H2,4'2, H^,^, and so on. So, by 
specifying Ho, 4>q and 4>i, all the remaining terms in the solution are fixed. The 
aim now is to choose the three input functions to ensure that successive terms 
in the series get progressively smaller. This provides just the right number of 
equations to specify all coefficients, save for two arbitrary coefficients in <pQ. 
This results in a series expansion controlled entirely by two arbitrary constants, 
which is the expected number of degrees of freedom once we have fixed the 
singularity to t = 0. In order to generate curvature it turns out that the input 
functions need to be power series in w 1 / 3 , which ensures that the scale factor 
goes as u 1 / 3 at small times. The series solution is only required to find suitable 
initial conditions for numerical evolution, so only the first few terms are required. 
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Expanding up to order u 5 / 3 we find that 




47r 2 2 4\/37r 2, \ 66567r6| 5/3 

3?r& | u 8/3 



1296tt 
and 



54V3^A - 216V3 7r 3/ V 2 &o + 36tt/x 2 6o) u 2 (26) 



1 



-V37r + 367r&o u. (27) 



12tt 216tt 

Under scaling the three key parameters in the model transform as fi h- > rr^, 
&o l—> &0j and ^4 l— * cr 4 / 3 &4. This scaling transformation for is entirely as 
expected, given that it is the coefficient of it 4 / 3 in the series for cf>Q. It follows 
that the quantity 64/ fi 4 ^ 3 is scale invariant. 

The fact that 64 controls the curvature can be seen from equation £Q| which, 
to leading order, yields 



5 = i(^! V'Y-^V' 2 , ■ ■ 

l p fi V125447T, 



-it &•»>■" + ■■•■ <»> 



Clearly, the arbitrary constant 64 must be negative for a closed universe. The 
terms on the right-hand side of equation l|28|l are all scale invariant, apart from 
the first factor of In figure [5] we illustrate Hubble function entry into 

the inflationary regime for typical values of 60 and 64//Z 4 / 3 of interest. We 
see that the onset of inflation corresponds to fiu ~ 0.1. In order to generate 
perturbations consistent with observation the scalar field must have a mass of 
the order of 10 _6 m p . It follows that the onset of inflation occurs at a time of 
around 10 5 Planck times. The radius of the universe at the onset of inflation is 
then given approximately by 

2 

R w — l p , (29) 

and so is of order 10 5 Planck lengths. Inflation therefore starts at an epoch 
well into the classical regime. Quantum gravity effects would be expected to 
be relevant when the radius of the universe is of the order of the Planck scale, 
which occurs when u ss [i 2 and is well before any inflationary period (for physical 
values of /Lt). We are therefore quite justified in running the evolution equations 
back past the inflationary regime, and right up to near the initial singularity. 
It is only when R = l p that the equations will break down, and we would 
look to quantum gravity to explain the formation of the initial, Planck-scale 
sized universe. Indeed, we would argue more strongly that we have to run the 
equations backwards in time to well before the start of inflation before reaching 
an epoch where new physics could be expected to enter the problem. 

As the universe is described by a 3-sphere of radius i?, the total volume of 
the universe is given to leading order by 

^ = 2 * 2 (il5;f <4^° + '--- (30) 
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Figure 9: Hubble function entry into and exit from the inflationary regime. 
The Hubble function emerges from the big bang going as l/3t. As H falls 
it quickly enters the inflationary region, during which it falls linearly. As the 
inflationary era ends the density and pressure start oscillating around values for 
a matter-dominated cosmology. The effective singularity for this dust cosmology 
is displaced from t = 0. The broken line plots the natural log of 2/3(t — 14.56). 
Time t is measured in units of t p , and H in units of t~ . The input parameters 
were set to bo = 2.48, 6 4 = —0.51 and m — m p (n= 1). 
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Figure 10: Evolution of the conformal time rj(t) as a function of\n{t). The 
parameters for the model are as given in figure |5] and again t is measured in 
units of the Planck time. This plot sets fj, = 1, so the time variable must be 
scaled to correspond to more physical values of /i « 10~ 6 . 



Following from this, an interesting calculation we can perform in a closed uni- 
verse is to find the total energy contained within it in the scalar field. By 
integrating the energy density we find that 

7T / 2187 \ 3/4 i h /~^/ 3 \ 3/2 h , , 

and when the scale-invariant quantity 64 //x 4 / 3 is or order unity (as required for 
physical models) we have 

Etot ~7^T- (32) 

This tell us that the action Et is very large, when measured in units of h. This 
is reassuring, as it means we are justified in treating the universe on the whole 
as a classical object. 

The plots in figure |5| illustrates the general behaviour of the Hubble function 
H . As the universe emerges from the big bang the energy density in the scalar 
field is dominated by the <f> term, and the field behaves as if it is massless. It 
follows that H initially falls as l/(3i). But once H has fallen sufficiently far we 
enter a region in which m 2 <fi 2 starts to dominate over cp 2 . These are suitable 
initial conditions for the universe to enter an inflationary phase. By varying 60, 
&4 and fi we control both the values of the fields as we enter the inflationary 
period, and how long the inflationary period lasts. The cosmological constant 
plays no significant role in this part of the evolution. The dynamics displayed 
in this figure is quite robust over a range of input parameters. A significant 
point here is that we enter the inflationary regime from a region of high H, as 
opposed to the H « value favoured in some models of quantum cosmology. We 
therefore never enter the regions of parameter space where the chaotic evolution 
noticed by Page [21] an d Cornish & Shellard [2H] is relevant. 

The typical behaviour as we exit the inflationary region is also shown in 
figure |5J The end of inflation is characterised by H ss [i. Beyond this point the 
scalar field enters an oscillatory phase, with the time-averaged fields satisfying 
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the conditions for a simple matter-dominated cosmology. In this case H can 
be approximated by a curve going as 2/3(< — to), describing a dust model with 
a displaced origin. The universe then appears as if it has been generated by 
an 'effective' big bang at a later time. Around this time we expect reheating 
effects to start to dominate, so that in reality the universe must pass over to 
a radiation-dominated era. But the naive 'effective big bang' concept is useful 
for extracting some qualitative predictions from our model @]. Figure ITUI shows 
the evolution of conformal time for the model depicted in figures ED As the 
universe emerges from the initial singularity, 77 grows at t 2 / 3 . But once the 
inflationary region is entered, R(t) starts to increase rapidly. So the conformal 
time, which involves the time integral of R , quickly saturates. So if rj has not 
reached ir/2 before R has inflated significantly the universe will have to exist 
for an extremely long time to reach the boundary value of n/2. Similarly, in 
this simple picture, we can obtain the prediction for A given in equation J5J) (for 
details see @j). We note with the number of e-folds N » 46, then A ~ 10 _122 Z~j 2 , 
which is in exactly the right ball park! 

5 The power spectrum 

In 0] we give an account of how to extract the power spectrum of primordial 
fluctuations, using the various approximations that have become standard in 
the literature. As an example of the type of results achieved, we show in figure 
HTI a comparison of our computed scalar power spectrum with the WMAP best 
fit power law and running spectral index fits. The vertical normalisations for 
these fits have been chosen arbitrarily, since it is the shape of the spectrum at 
large k which is of greatest interest here. The running spectral index graph is 
interesting in that it suggests that this model is attempting to emulate both the 
cutoff at low k, as well as the reduced power at large k, which occurs in our 
model. The latter occurs as a result of computing the perturbation spectrum 
more accurately than usual in this regime, and is not a specific feature of using 
a closed universe model. 

In figure IT21 we show a predicted CMB power spectrum from our model, in 
comparison with the WMAP results and the prediction from a strictly power 
law initial spectrum. This model has f^o = 1-04, so is consistent with WMAP 
at the la level. We also take D, h h 2 = 0.0224, h = 0.60 and ft cdm /i 2 = 0.110. 
Together these yield a value of f2 c dm = 0.306, which is reasonable, although 
of course the Ho value is rather low. It can be seen that our predicted CMB 
spectrum is in much better agreement with the WMAP results at low I than 
the predictions from the primordial power law. This could be taken as good 
evidence for our model, except that it is necessary to check the correctness of the 
approximations we have used in calculating the perturbation power spectrum. 
This work is currently in progress, but we end this paper with a brief discussion 
of some of the technical problems involved. 

The starting point for calculating the scalar curvature spectrum in a simple 
flat model is to write the linearised perturbed action in the form 

5 2 S = i J drj^x [v a - ti^v^vj + % 2 j . (33) 

Here v is a gauge-invariant combination of matter and metric perturbations, 
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Figure 11: Comparison of the scalar power spectrum of curvature perturbations, 
V-n(k), with power law models. The function plotted is lO w (2ir) 2 V-ji(k) as a 
function of log 10 (fc), assuming h = 0.65. The solid line represents the numerical 
predictions from our model. The long dashes represent the best fit power law 
(n s = 0.96) and the short dashes are the WMAP running spectral index best 
fit. Notice that the vertical scale runs from 800 to 1300, so the differences are 
slightly exaggerated. 




I 

Figure 12: CMB power spectrum for a model with Qq = 1.04 The parameters 
are discussed in the text. The experimental points shown with la error bars are 
the WMAP results £Q and the dashed curve corresponds to the best fit A CDM 
power law CMB power spectrum as distributed in the WMAP data products. 
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dashes denote the derivative with respect to conformal time, and z is defined by 
z = (j)' /H, where 4>o is the unperturbed field. By working with this action the 
entire problem is reduced to one of analysing a scalar field with a time-dependent 
mass term. This approach generalises to the non-flat case, although here the 
definition of z becomes more complicated [22] ■ Also, for closed models, the 
mode expansions necessary to carry out quantisation have to be performed using 
spatial sections given by the 3-sphere S 3 , so that the 'comoving wavenumber' 
k takes on integer values, with k = 3 its lowest (non-gauge) value. The mode 
equations found this way are very complicated, but retain the general form 

v'k + {k 2 ~ f(v, k)) v k = 0, (34) 

where in the flat case f(rj, k) would be z"/z, and in all cases f(rj, k) is calculable 
from knowledge of the background evolution. 

The challenge now is to find suitable 'quantum initial conditions' so that after 
evolution through inflation, the variables v k can be used to find the perturbation 
spectrum. To achieve this involves a mode decomposition of v k into positive and 
negative frequencies. The standard way to approach this is to assume that in 
the asymptotic past the background is either Minkowski or de Sitter, so that 
one knows the correct vacuum to chose. But this is clearly inappropriate here, 
as looking back in time we encounter the initial singularity. 

The question of how to proceed in the absence of any asymptotic notion of 
the vacuum state has been widely discussed. One simple approach is provided by 
Hamiltonian diagonalisation, where on each time slice modes are decomposed 
into positive and negative energy states by the Hamiltonian operator. But 
this technique tends to overestimate the particle production rate. A clear way 
to proceed was developed by Parker and Fulling j^U EH] and introduces the 
concept of an adiabatic vacuum. The application of this to the present case 
is complicated, but feasible, and initial results support a low-Zc cutoff, though 
with less pronounced effects than we have found using the standard approximate 
techniques. This will be the subject of a future publication. 
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